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We derive shape-independent limits to the spectral radiative heat-transfer rate between two closely spaced 
bodies, generalizing the concept of a black body to the case of near-field energy transfer. Through conservation 
of energy and reciprocity, we show that each body of susceptibility x can emit and absorb radiation at enhanced 
rates bounded by |x//Inix» optimally mediated by near-field photon transfer proportional to 1/d^ across a 
separation distance d. Dipole-dipole and dipole-plate structures approach restricted versions of the limit, but 
common large-area structures do not exhibit the material enhancement factor and thus fall short of the general 
limit. By contrast, we find that particle arrays interacting in an idealized Born approximation (i.e., neglecting 
multiple scattering) exhibit both enhancement factors, suggesting the possibility of orders-of-magnitude im¬ 
provement beyond previous designs and the potential for radiative heat transfer to be comparable to conductive 
heat transfer through air at room temperature, and significantly greater at higher temperatures. 


Heat exchange mediated by photons, or radiative heat trans¬ 
fer, can be dramatically modified for bodies separated by 
small gaps [1-7]. We exploit energy-conservation and reci¬ 
procity principles to derive fundamental limits to the near¬ 
field spectral heat flux between closely spaced bodies of ar¬ 
bitrary shape, given only their material susceptibilities x(<^) 
and their separation distance d. Our approach enables us to 
define optimal absorbers and emitters in the near field, which 
contrast sharply with far-field black bodies: their response is 
bounded by the amplitude of their volume polarization cur¬ 
rents, rather than their surface absorptivities, and maximum 
energy transfer requires coordinated design of the two bod¬ 
ies (whereas the far-field limit derives from the properties 
of a single black body). These distinguishing characteris¬ 
tics lead to two possible enhancements relative to black-body 
emission: a material enhancement factor |x(cc)p/Imx(cc) 
that represents the maximum absorber and emitter polar¬ 
ization currents, and a near-field enhancement factor 1/d^ 
that represents maximum interaction between currents in free 
space. We show that restricted versions of our limits can 
be approached for sphere-sphere and sphere-plate configu¬ 
rations. For two extended structures, however, common pla¬ 
nar geometries—including bulk metals [8-17], metamateri¬ 
als [18-24], and thin films [25-31] —exhibit flux rates orders 
of magnitude short of the limits because they do not satisfy 
the optimal-absorber condition. Instead, we find that idealized 
plasmonic-particle arrays, interacting within a Born approxi¬ 
mation with negligible multiple scattering, approach the limits 
at selected frequencies, and that the possibility of reaching the 
limits, even over a narrow bandwidth (a desirable feature for 
thermophotovoltaics [7, 32-35]), would represent an orders- 
of-magnitude improvement over current designs. 

A ray-optical black body absorbs every photon incident 
upon its surface, which by reciprocity (Kirchoff’s Law) yields 
its emissivity and the black-body limit to thermal radia¬ 
tion [36]. At wavelength and subwavelength scales, nanos¬ 
tructures can exhibit optical cross-sections much larger than 
their physical cross-sections [37], making it difficult even to 
define quantities like emissivity. A further difficulty in the 
near field is the presence of evanescent waves, which can in¬ 


crease transmitted power but only through interference with 
reflected waves [38]. Although the possibility of enhancement 
beyond the blackbody limit was realized by Rytov, Polder, and 
others in the 1950s [1,2], efforts to find underlying limits have 
been restricted to planar structures with translation symmetry 
(including metamaterials), without consideration of material 
loss [10, 15-17, 21, 30]. Spherical-harmonic [39, 40] and 
Green’s-function [41] limits are difficult to apply in the near 
field where a large but unknown number of spherical harmon¬ 
ics can be excited by general shapes [42] . 

Without reference to particular structures or symmetries, 
assuming only linear electromagnetism, we translate the reci¬ 
procity principle to the near field by applying it to polarization 
currents within the bodies. Dipoles in vacuum exchange en¬ 
ergy at a rate limited by the energy density of an outgoing free- 
space wave [43]. As we show below, the maximum energy 
transfer between material bodies occurs when the currents 
within the bodies couple individually at the dipole-dipole 
limit, amplified by material enhancement factors. These con¬ 
ditions allow for much greater heat transfer than has previ¬ 
ously been shown possible. 

Radiative heat exchange is depicted schematically in 
Fig. 1(a): fluctuating currents arise in body 1 at temperature 
Ti, and transfer energy to body 2 at a rate of [4] 

POO 

= / Hoj) [0(cu, Ti) - 0(u;, T 2 )] dec, (1) 
Jo 

where <h(cc) is a temperature-independent energy flux and 0 
is the Planck spectrum. <F(cc) is the designable quantity of 
interest, to be tailored as a function of frequency depending 
on the application and available materials. 

Limits —The spectral heat flux 4>(cc) is the power absorbed 
in body 2 from fluctuating sources in body 1 (or vice versa). In 
recent work [42] we have bounded the scattering properties of 
any dissipative medium excited by a known, externally gener¬ 
ated incident field. The bounds arise from the functional de¬ 
pendencies of power expressions with respect to induced cur¬ 
rents: absorption is a quadratic functional, whereas extinction 
(absorption-Fscattering), given by the optical theorem [44-47], 
is only a linear functional. Energy conservation requires that 



2 



FIG. 1. (a) Radiative heat transfer: Fluctuating currents in an emitter (body 1, susceptibility xi) generate a field Einc,i and transfer energy to 
an absorber (body 2, susceptibility X 2 ) at a rate Pabs, 2 . (b) Energy conservation bounds Pabs ,2 in terms of Einc,i, and a resonant enhancement 
factor 1x21^ / Im X 2 . (c) The sources and “receivers” can be exchanged by reciprocity, whereupon (d) absorption in body 1 is bounded, yielding 
a spectral-fiux limit determined by Xi’ X 2 , and the free-space GF Gq. For near-field transfer the GF integral is ~ 1/d^, for separation d. 


extinction be greater than absorption, which imposes a bound 
on the magnitude of the excited currents. Radiative heat trans¬ 
fer, however, involves sources within one of the scatterers, 
preventing a simple optical theorem. 

To circumvent this issue we reframe the scattering prob¬ 
lem (without approximation). We define the “incident” field 
to be the unknown field emanating from body 1, and the 
“scattered” field to arise only with the introduction of body 
2. For a Green’s function (GF) that is the field of dipole 
in the presence of body 1, the fields are given by a standard 
integral-equation separation [48], Einc,i = {i/sQUj) fy^ GiJ 
and Escat,i = fy^ GiP, where J are the stochastic source cur¬ 
rents in body 1, P is the polarization field induced in body 2, 
and £0 is the vacuum permittivity. This decomposition permits 
an optimal theorem with respect to body 2, such that its extinc¬ 
tion is proportional to Im fy^ Einc,i • P (its absorption [44] is 

proportional to fy^ |P|^)- The energy-conservation arguments 
from above imply that absorption in body 2 is bounded, 


P abs ,2 ^ 


epu) |X2(^)P 
2 Iinx2(w) 



( 2 ) 


which is formally derived by variational calculus [42]. To 
achieve this limit, the optimal polarization field must be pro¬ 
portional to the incident field, P ^ Einc,i, to maximize the 
extinction overlap integral. In the near field, where source 
fields rapidly decay, negative-permittivity metals that support 
surface-plasmon modes can achieve this condition, as we will 
demonstrate. 

The limit in Eq. (2) reduces the optimal-flux prob¬ 
lem to a question of how large the emitted fleld 
Einc,l can be in V 2 . Inserting Einc,i into Eq. (2) 
yields an integral of the stochastic currents, which is 
determined by the fluctuation-dissipation theorem [4], 

(jj(x,w), Jfe(x',w)) = 4eoW0(w,Ti)Iin[x{a;)](5jfe(5(x - 
x')/7r, such that the ensemble-averaged emitted fleld at X 2 in 
V 2 is (|Ei„c,i(x 2 )|^) = 4eow0(Imxi)/i/^ ||Gi(x 2 ,xi)||^, 
where H-H^ denotes the Erobenius norm [49]. By reci¬ 
procity [50] one can exchange the positions in the integrand, 


xi ^ X 2 (while transposing the GE, but the transpose does 
not affect the norm), such that emission from Vi is equivalent 
to absorption for free-space sources in V 2 , as in Eig. 1(c). Ab¬ 
sorption is bounded by energy conservation [42], limiting the 
emitted-fleld magnitude: 


|Ei, 


:,1(X2)|^) < 4£oW0i^^^ / ||Go(xi,X2)||^ (3) 

' ImXi Jxi 


where Gq is free-space GE, cf. Eig. 1(d). Inserting Eq. (3) 
into Eq. (2) and separating the Planck spectrum by Eq. (1), the 
maximum flux between two bodies is 


T>(u;) < 


2 IxiHI^ \X 2 Hf f f 

TT Iinxi(w) Iinx2(w) Xj iv2 


|Go(xi,X2) 


IF • 


( 4 ) 


The limit of Eq. (4) can be further simplifled. In the near fleld, 
Go is ideally dominated by the quasistatic term ^ 1/r^, which 
is primarily responsible for the evanescent waves that en¬ 
able greater-than-black-body heat-transfer rates [4, 7]. Drop¬ 
ping higher-order terms (further discussed in [51]), we bound 
Eq. (4) by integrating over the inflnite half-spaces contain¬ 
ing Vi and V 2 , assuming a separating plane between the two 
bodies. (If not, e.g. between two curved surfaces, only 
the coefficients change.) Eor bodies separated by a distance 
d, the integral over the (infinite) area A is given by [51] 
•fvi,V2 |||. = A/32F(i^, yielding flux limits per area or 
relative to a black body with flux T>bb = A:^A/4f^ [4]: 

< 1 lxi(^)l^ |X 2 (^)|^ 

A ~ 167r^(i^ Imxi(<^) Imx2(<^) 

^ 1 IxiHI^ Ix2HI^ 

^bb(w) “ 4(A:d)2 Imxi(w) Imx2(w)' 

Eqs. (4-6) are fundamental limits to the near-field spec¬ 
tral heat flux between two bodies and form the central re¬ 
sults of this Letter. They arise from basic limitations to the 
currents that can be excited in dissipative media, and their 
derivations further suggest physical characteristics of the op¬ 
timal response in near-field heat transfer: an optimal emitter 
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enhances and absorbs near-field waves from reciprocal exter¬ 
nal sources in the absence of the absorber whereas an opti¬ 
mal absorber enhances and absorbs near-field waves from the 
emitter, in the presence of the emitter. These principles can 
be understood by working backwards through Fig. 1. The 
optimal-emitter condition identifies the largest field that can 
be generated in an exterior volume (V 2 ) by considering the 
reciprocal absorption problem, per Fig. 1(c). Reinserting the 
absorber, cf. Fig. 1(b), should not reflect the emitted field 
but rather enhance and absorb it. Because heat flux is sym¬ 
metric with respect to absorber-emitter exchange, both bod¬ 
ies should satisfy each condition (induced currents propor¬ 
tional to source fields). Eq. (4) can be interpreted as sources 
throughout the emitter generating free-space dipolar fields 
enhanced by |xi|^/Imxi, which are further enhanced by 
|X 2 |^/IniX 2 and absorbed. The dipole-dipole interactions 
are bounded by their separation distance [43, 52], leading to 
simple shape-independent limits in Eqs. (4-6). Ideal struc¬ 
tures that achieve these limits can have significantly greater 
heat transfer than black bodies, even if their spectral fiux has 
a narrow bandwidth. Whereas the heat transfer between black 
bodies in the far field H/A = where ctsb is the 

Stefan-Boltzmann constant [36], a straightforward calcula¬ 
tion citeNotel shows that ideal near-field heat exchange over 
a narrow bandwidth Alo/uj = Imx/lxl^ typical of plasmonic 
systems [53, 54], can achieve per-area transfer rates of 


2 ixi^ 

A 7{kd)‘^lmx' 


(7) 


exhibiting both distance and material enhancements relative 
to the Stefan-Boltzmann rate. 

The limits generalize [51] to local media with ten¬ 
sor susceptitbilities via the replacement |xP/Imx ^ 

x(Imx)~^X^ . Nonlocal effects, which appear below 

2 

lOnm length scales [55] and which regularize the l/cP diver¬ 
gence [4], are outside the scope of these limits, but we believe 
that a generalization to nonlocal x is possible and have pre¬ 
liminary results [56] suggesting that “hydrodynamic” [57, 58] 
nonlocal materials cannot not surpass the local-x bounds. 

Dipolar Interactions —If one of the bodies is small enough 
for its response to be dipolar, the optimal-absorber and 
optimal-emitter conditions converge: the polarization currents 
induced in each structure by free-space dipoles in place of the 
opposite structure must be proportional to the incident fields. 
This condition is satisfied for two-dipole transfer, and the en¬ 
hancement of the emitted and absorbed fields is possible via 
“plasmonic” resonances in metallic nanoparticles. For two 
identical particles with volumes V, tip-to-center-of-mass dis¬ 
tances r, and tip-to-tip separation d, Eq. (4) limits the fiux: 


[^(^)] dipole-dipole — 


3 |x2MI^ 

47r3 Iinxi(w) Iinx 2 {w) (2r + df 


( 8 ) 


The radiative fiux between quasistatic metal spheres is known 
analytically [4] and peaks at the limit given by Eq. (8). 



FIG. 2. Comparison of heat flux in sphere-sphere and sphere-plate 
structures to the analytical limits of Eqs. (8,9). Two Drude metal 
spheres (orange circles, fit to a solid line) approach the dipole-dipole 
limit (dashed orange) at their resonant frequency, Ures ~ ojp /V3. A 
sphere and a plate (blue circles) approach within a factor of two of 
the limit between dipolar and extended objects (dashed blue), if the 
material resonance of the plate is slightly modifled (see text). In each 
case the separation is d = 0.1c/Ures, with sphere radii r = d/5. The 
flux rates exhibit the material enhancement factor |x|^/x)^^ but 
not the near-fleld enhancement factor, due to the lack of large-area 
interactions. The sphere area A is taken to be the cross-section nr^. 


Heat transfer between a dipole and an extended structure 
is limited by integrating over the half-space occupied by any 
extended structure, yielding a maximum flux 


[^(^)]dipole-to-ext — 


1 lxi(^)l^ Ix2HI^ V 

Stt^ Imxi(<^) Imx2(<^) + d)^ ’ 


(9) 


where r + d is the distance between the extended structure and 
the particle’s center. Heat flux between a sphere and a bulk 
metal, each supporting a plasmonic mode, can achieve half of 
the maximum flux [4, 51, 59] if the resonances align. This 
geometry falls short by a factor of two because planar surface 
plasmons exist only for TM polarization [60], and thus the pla¬ 
nar structure reflects near-field TE-polarized light emitted by 
the sphere. Neither structure exhibits the 1/d^ enhancement 
factor, which for dipolar coupling (^ 1/d^) requires interac¬ 
tions over two extended areas. 

Fig. 2 compares flux rates for sphere-sphere (orange cir¬ 
cles) and sphere-plate (blue circles) geometries, computed 
by the fluctuating-surface current method [61-63], to the 
limits of Eqs. (8,9) (orange and blue dashed lines, resp.). 
The spheres are modeled by Drude susceptibilities [44] with 
plasma frequency uOp and loss rate 7 = O.lujp. The “plate” is 
simulated by a very large ellipsoid (volume ^ 7000 x larger 
than the sphere) comprising a material with a modified plasma 
frequency, = ^^^JSujp, and a modified loss rate, 7 pi = 
27/3, to align the resonant frequencies of the sphere and plate 
without modifying the flux limit. In each case the separation 
distance d = O.lc/ccres and the sphere radii are r = d/b. 
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FIG. 3. (a,b) Comparison of heat flux between mirror images of large-area Drude-metal structures separated by d = 0.1c/ct;p. (a) Structures 
optimized for maximum flux at three frequencies, u = (0.2,0.4, l/\/2)cjp, for a material loss rate 7 = O.Olcjp. Thin Aims (purple), 
hyperbolic metamaterials (blue), and elliptical metamaterials (orange) exceed black-body enhancements but fall far short of the limit (black) 
from Eq. (5). The dashed silver line represents the heat transfer for an idealized plasmonic-particle array without multiple scattering, (b) 
Optimized structures as a function of loss rate, for u = 0.4cjp. Each structure exhibits the 1/d^ near-fleld enhancement factor, but only the 
idealized particle array exhibits the \x \^/(to x)^ ~ 1/7^ material enhancement factor, (c) Erequency-integrated heat transfer coefficient of a 
structure that reaches the single-frequency limit in Eq. (5) over a narrow bandwidth Auj oc 7 . Radiative heat exchange in this limit shows the 
possibility of surpassing conductive heat transfer through air (dotted) at T = 300K (gold), which is not possible for plate-plate conflgurations 
(inset, dashed), and of signiflcant further enhancements at higher temperatures (blue, purple). 


The computations support the analytical result that the dipolar 
limits can be approached to within at least a factor of two. 

Extended Structures —For extended structures that do not 
behave like single dipoles, the optimal-absorber constraint is 
more demanding in that the absorber should enhance the emit¬ 
ted field while accounting for interactions between the two 
bodies. We will show that common planar structures do not 
exhibit this behavior but that nano structured media offer the 
possibility of approaching it. 

Bulk metals (negative-permittivity materials) support sur¬ 
face plasmons that enable greater-than-blackbody heat flux at 
their resonant frequency. Individually, a single metal inter¬ 
face nearly satisfies the optimal-emitter condition, emitting 
near-field waves over a broad bandwidth of surface-parallel 
wavevectors (which enabled the nearly optimal sphere-plate 
transfer above). However, when a second metal is brought 
close to the first, it reflects most of the incident field, ex¬ 
cept over a narrow wavevector-bandwidth, due to multiple¬ 
scattering effects between the bodies. The failure of the 
two-metal geometry to achieve the optimal-absorber condi¬ 
tion leads to a peak spectral heat flux, at the surface-plasmon 
frequency ccgp, of approximately citeNotel 


^(^sp) 

- ^ In 

Ixl" 1 

A 

47r^d^ 

plate-to-plate ' ^ 

_4(Imx)^_ 


which is significantly smaller than the limit in Eq. (5) due to 
the weak, logarithmic material enhancement. Eq. (10) appears 
to be new and is a significantly better approximation than pla¬ 
nar bounds that do not account for material loss [10, 16], as 
discussed in the SM citeNotel. The shortcomings of the bulk- 
metal interactions cannot be overcome with simple metamate¬ 
rial or thin-film geometries. The fiux rate between hyperbolic 
metamaterials (HMMs) is material-independent [21, 51]. Op¬ 
timal thin films behave similarly to HMMs [31], thereby also 


falling short of the limits. “Elliptical” metamaterials, with 
nearly isotropic effective permittivities, exhibit resonances for 
Xeff ~ — 2 and thus transfer heat at a rate similar to Eq. (10), 
limited by the same interference effects discussed above, and 
because |xeff|^ « Ixl^- 

Eig. 3(a,b) demonstrates the shortcomings of such struc¬ 
tures, showing the computed heat fiux between mirror im¬ 
ages of thin-film (purple), hyperbolic-metamaterial (blue), 
and elliptical-metamaterial (orange) structures, as a function 
of (a) frequency and (b) material-loss rate, for a fixed sep¬ 
aration d = O.lc/ccp. Assuming smooth surfaces without 
roughness, the structural parameters are computationally op¬ 
timized citeNotel using a derivative-free local optimization 
algorithm [64, 65]. Eig. 3(b) shows that the sub-optimal 
performance can be attributed primarily to the fact that the 
structures do not exhibit the material enhancement factor 
1x7 (Imx)^ 1 / 7 ^, as predicted by Eq. (10) and due to 
the significant reflections in such geometries. 

The spectral heat fiux of the limit in Eq. (4) can be inter¬ 
preted as the exchange of enhanced free-space dipole fields, 
as discussed above. Guided by this intuition, we include 
in Eig. 3(a,b) the heat fiux between close-packed arrays of 
oblate disk ellipsoids (dashed silver lines), small enough to 
be dipolar. We idealize their response as the additive sum of 
Eq. (8) over a lattice neglecting multiple scattering (i.e. in a 
Born approximation) [66] and accounting for the polarization- 
dependence of non-spherical ellipsoids [37]. This structure 
combines the individual-particle interactions that exhibit the 
material enhancement (which planar bodies do not) with the 
large-area interactions that exhibit l/d‘^ near-field enhance¬ 
ment (which isolated bodies do not). Eig. 3(a,b) suggest the 
possibility for two to three orders of magnitude enhancement 
by periodic structuring and tailored local interactions. 

Experimental measurements of radiative heat transfer are 
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done in vacuum [ 6 , 12 , 13] because radiative transfer is domi¬ 
nated by conductive transfer through an air gap. Achieving the 
limits presented here, even over a narrow bandwidth, could 
transform this landscape. Fig. 3(c) shows the heat-transfer 
coefficient h = f ^(dO/dT)dcj for extended Drude-metal 
structure with loss rates 7 = O.Olujp (appropriate e.g. for 
Ag and Au [67]). For Lorentzian-shaped energy transfer with 
tunable center frequency ujj-es = peaked at the limit 

given by Eq. (5), with a bandwidth Acc = 7 [51, 53, 54], ra¬ 
diative transfer can surpass conductive (thermal conductivity 
i^air = 0.026W/m-K [ 68 ]) even at T = 300K. In the inset 
we fix the wavelengths at A = 7.6/im for T = 300K and 
A = 3/im for T = ISOOAT, and plot as a function of distance 
for plate-plate (dashed) and optimal (solid) transfer. We find 
that radiative transfer can surpass conductive at separation of 
d = 50nm at 300K and almost d = 0.5/im at T = 1500K, 
gap sizes that are readily achievable in experiments. 

Radiative heat transfer at the nanoscale is a nascent but 
growing field. Calculations have primarily been for dipo¬ 
lar [5, 10, 59] or highly symmetric bodies [8-14, 18-23, 25- 
31, 69-71], with computational study of more complex ge¬ 
ometries possible only recently [35, 62, 63, 72-74]. We have 
show that, guided by the physical principles presented here, a 
targeted search through the mostly uncharted near-field design 
space offers the prospect of orders-of-magnitude enhance¬ 
ments in radiative energy transfer. 
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We provide: (1) a derivation and discussion of the higher-order terms in the heat flux limits, which tend to 
be very small for near-fleld heat transfer, (2) a derivation of the peak heat flux between two Aims, for a given 
material, (3) a derivation of the limits for a very general class of materials, (4) a derivation of the radiative heat 
transfer coefficient if the limiting flux rates are achieved, alongside a comparison to conductive heat transfer 
through air, and (5) an estimate of the frequency-integrated heat transfer for a narrow-band resonance. 


I. EVALUATION OF INTEGRAL LIMITS AND 
HIGHER-ORDER TERMS 

In this section, we present calculations and clarify the step 
needed to go from Eq. (5) to Eq. (6) of the main text. Specif¬ 
ically, Eq. (5) is an integral bound that applies to any near- or 
far-held interactions, depending only on conservation of en¬ 
ergy arguments. Eq. (6) simplifles the bound for the case of 
near-fleld heat transfer by assuming that the near-field qua¬ 
sistatic 1/r^ term in Gq is the dominant term and integrat¬ 
ing over the infinite half-spaces occupied by the two bodies. 
(All equations and figures in this Supplementary Material are 
preceded with an “S,” whereas equations and figures with¬ 
out an “S” refer to the main text.) Here we justify dropping 
the 1/r^ and 1/r terms in the Green’s function. Although 
for many structures it is known that optimal near-field heat 
transfer is governed by high-wavevector waves corresponding 
to the 1/r^ term, the mathematical justification for dropping 
the terms is somewhat subtle. Integrated over infinite half¬ 
spaces, the two terms diverge. We show that this divergence 
is unphysical—originating from the optimal variational fields 
that are appropriate in the near field but which do not satisfy 
Maxwell’s equations in the far field. Moreover, we show that 
for finite, reasonable interaction distances, their contributions 
are negligible compared to the contribution of the 1/r^ term. 
As shown in the text, the limit of Eq. (5), keeping only the 
1/^3 ^erm, yields very good agreement with the response of 
sphere-sphere and sphere-plate interactions. 

The squared Erobenius norm of the homogeneous Green’s 
function is: 


l|Go||^ = 


Stt^ 


3 1 1 

{kr)^ {kr)^ {kr)‘^ 


(S.l) 


which has contributions from 1/r^, 1/r^, and 1/r^ terms. Eor 
convenience, instead of taking infinite half-spaces, we assume 
that both bodies are contained within a circular cylinder of 
radius R and height L. The integral of the norm over both 
volumes is a six-dimensional integral, but we bound it above 
by fixing the source in one body at its center (x = y = 0), and 
multiplying by the cylindrical area A = : 



J dzi j dz2 j dp2'Kp\\Gi^\ 


(S.2) 


finitely wide) structures with translational and rotational sym¬ 
metry; since we are interested in global bounds encompassing 
large structures it is thus a good approximation. The bound in 
Eq. (6) of the main text comes from the 1 /r^ term in the GE 
for an infinite volume (it is very weakly decreased for large 
but finite structures). The integral is given by: 
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where V-[ and V 2 are the infinite half-spaces containing the 
bodies. Multiplying by the prefactors in Eq. (S.l) yields the 
bound in Eq. (6) of the main text. Over finite volumes, the 
second term is more complicated: 
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The third term is given by: 
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Eqs. (S.4,S.5) are difficult to disentangle so we consider large 
but finite volumes. Large bodies satisfy 

L,R'>d (S.6) 


where we have further simplified the integral using cylindri- such that their sizes are much larger than their spacing. Not 
cal coordinates. The multiplication by A is exact for (in- only do L and R represent the physical sizes of the bodies, 

































2 


they also represent the interaction sizes: they are the volumes 
over which polarization currents within the respective bod¬ 
ies transfer energy. Near-field interactions by definition oc¬ 
cur between charges or currents at the subwavelength scale, 
such that one is typically interested in sizes L X. Con¬ 
versely, surface waves between structures are example of co¬ 
herent subwavelength interactions that potentially take place 
over distances much greater than the wavelength, i? ^ A. 
Thus the finite-but-large asymptotic expansion relevent for 
near-field heat transfer can be made by taking 

d < L < i? (S.7) 

for the two circular cylinders with radii R, heights L, and sep¬ 
aration distance d. In this asymptotic limit, the terms simplify: 

jy.y. ^ = 

The divergences in the second and third terms are relatively 
weak. The second term is negligible compared to the third 
term, which tends to be very small compared to the first. 
The comparison between the first and third term essentially 
compares l/{kd)‘^ versus {kL)‘^; even in a generous upper 
bound in which kL ^ 1, the third term is still much smaller 
than l/{kd)‘^ ^ 1. In Table 1 we compare the bound aris¬ 
ing from Eq. (6) to the bound that would arise from adding 
Eqs. (S.4,S.5) to Eq. (6). We see that for near-field distances 
(d X), even very large estimates of the interaction distances 
L and R lead to only small modifications to the upper limit, 
on the order of 1% and in some cases significantly smaller. 


kd 

kL 

kR 

Eq. (6) 

Eq. (6)+Eqs. (S.4,S.5) 

Rel. Error 

0.01 

1 

1 

1250 

1252 

0.17% 

0.01 

1 

10 

1250 

1254 

0.35% 

0.01 

1 

100 

1250 

1256 

0.53% 

0.001 

1 

100 

1.25 X 10® 

1.25008 X 10® 

0.0063% 

0.001 

10 

1000 

1.25 X 10® 

1.255 X 10® 

0.38% 


Einally, we note that these divergences arise even for far- 
field interactions, where they are clearly unphysical because 
finite blackbody limits to the flux per unit area are well known. 
The unphysical divergences arise from the assumption that 
the optimal polarization fields are proportional to the incident 
fields. Such a condition is ideal and achievable for the 1/r^ 
contribution of Go that typically dominates near-field trans¬ 
fer, but is unphysical for the more slowly decaying 1/r^ and 
1/r terms: a constant energy fiux is maintained in a lossy 
medium over large length scales, which is physically impossi¬ 
ble. One approach would be to “split” the problem into near- 
and far-field contributions, and to bound the interactions sep¬ 
arately. However, given the relatively weak nature of these 



FIG. SI. Heat flux per unit area of two Drude-metal bulk media 
as a function of material loss rate, at the resonant frequency 

uj — (jjply/2 and at a fixed separation of d = O.lc/ct;. Except 
for very large loss, the heat flux approaches the approximate rate 
of Eq. (S.22), confirming the logarithmic dependence on the mate¬ 
rial loss rate. Conversely, the limit of Ref. 2 is overly optimistic, and 
the “limit” of Ref. 3 is overly pessimistic. For the interatomic spac¬ 
ing that enters the limit of Ref. 2, we took b/X ^ 1/1000, which is 
appropriate e.g. for silver. 

contributions for finite interaction distances (< 1%), they can 
be ignored for near-field radiative heat transfer, justifying the 
use of Eq. (6) in the main text. 


II. HEAT TRANSFER BETWEEN BULK PLANAR MEDIA 

We derive the optimal heat-transfer rate between two pla¬ 
nar bodies comprising a material of susceptibility x(cc), cor¬ 
responding to Eq. (10) of the main text. Ref. 1 assumed a 
frequency-independent susceptibility, which they optimized 
for maximum heat transfer, whereas we assume a fixed (pos¬ 
sibly frequency-dependent) susceptibility. Ref. 2 and Ref. 3 
also provide expressions for optimal heat fiux between planar 
bodies, but their limits require wavevector-dependent mate¬ 
rial properties. The limits in both Ref. 2 and Ref. 3 arise only 
because a finite maximum surface-parallel wavevector magni¬ 
tude (/c||) is postulated: in Ref. 2 the maximum /cy ^ax = 1/^ 
is chosen, where b is the interatomic spacing of the metal; in 
Ref. 3, the maximum k\\ is inversely proportional to the gap 
spacing d, which does not account for large wavevectors that 
are possible when material losses are small. Although the in¬ 
teratomic spacing certainly sets an upper bound to the process 
as described by bulk materials, for lossy materials the loss is 
the limiting factor, not the interatomic spacing. We find a log¬ 
arithmic dependence (and divergence) of the heat fiux with the 
material loss rate, which we validate in Eig. SI. 

The radiative heat fiux 4>(u;) between two planar slabs is 
given by [4] 

A 

^ ^ ^11 


(S.ll) 
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where A is the area of the plates, is the magnitude of the 
surface-parallel part of the wavevector, and and Tp repre¬ 
sent the field transmissions from slab 1 to slab 2 for s and p 
polarizations, respectively. By symmetry, the surface-parallel 
wavevector is a conserved quantity between plane waves 
in each medium. The heat flux is characterized by a strong 
peak at a single k\\ (for a given cc) corresponding to the metal- 
insulator-metal plasmonic mode. We show that at a given fre¬ 
quency, the bandwidth in k\\ is approximately constant, while 
the peak energy transmission scales logarithmically with the 
inverse of the material loss rate. 

In the near field, we can focus only on the p-polarized trans¬ 
mission coeffient for evanescent waves with k\\ > ko. As¬ 
suming two slabs of the same material, with reflectivity r for 
waves incident from air, the transmission coefficient is [4] : 


4[Im(r)]^e 
|1 — 


(S.12) 


form at the moment, but we will note that for the optimal sus¬ 
ceptibility (see below) the real part of r is 0 and the imaginary 
part is potentially large. It we define the average [weighted by 
f{x)] value of x ns xq, it follows that f xf{x) = xq f f{x) 
and hence ^ can be approximately given by: 

The integral of / can be worked out: 



[Im(r)]^ TT /1 — Re(r^) \ 

Im(r^) 2 \ Im(r^) ) 

[Im(r)]^ _,f Im(r^) \ 
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(S.18) 


where 7 = k\\^l — /cg/Zcy ^ k\\, assuming k\\ ^ ko. With¬ 
out the denominator, Eq. (S.12) would yield a |x|^/(Imx)^ 
enhancement from the plasmon waves at each surface, man¬ 
ifested in the poles of Imr^ [5]. However, at the small 
distances necessary to transfer energy, the denominator— 
heuristically originating from the infinite sum of reflected 
waves—has an identical pole that cancels the one in the nu¬ 
merator. The resonances of Tp are instead metal-insulator- 
metal modes, with energy levels split around the single¬ 
surface plasmon energies [ 6 ], as discussed in the main text. 

Ref. 2 and Ref. 3 find limits to the transfer by noting that at 
every k\\ the maximum value of T is 1 (note that for conven¬ 
tional metals such a tranmission would require a wavevector- 
dependent permittivity). They define 
^p,max = [3], respectively, yielding limits: 

A 

■^(cc) 

A 

for interatomic spacing b and separation distance d. 

Instead we seek a limit assuming a conventional 
(wavevector-independent) material susceptibility x{^)- 
Defining x = 2k\\d, the flux is given by: 

A [Im(r)]^xe“^ 

Jo 1 — 2 Re(r^)e“^ + 

A 

where the integral lower bound can be set to zero because we 
have assumed kod 1 , and f{x) is defined by 

At large k\\, the reflectivity r is approximately constant and 
given by r = {e — l)/{e T 1). We will not insert its exact 


max,Ref. 2 


max.Ref. 
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(5.13) 

(5.14) 


where we used tan ^{l/x) = 7r/2 — tan ^(x), and for 
X small, tan“^(x) ^ x. For the final step, we can write 
Re(r^) = [Re(r)]^ — [Im(r)]^ = 1 — [Im(r)]^. To find the 
value of Xq, we approximate it (verifying later) as the value of 
X at which f{x) peaks. Setting the derivative of / in Eq. (S. 16) 
to zero yields: 


xo = In |rp. (S.19) 

Because r = 1/[1 + 2/x(cc;)], the optimal frequency for max¬ 
imum |r| is given by the frequency such that Re( —l/x(cc;)) = 
1 /2. At this frequency, r = i|xp/2Imx and we have: 



(5.20) 

(5.21) 


Thus at the optimal frequency, maximum energy transmission 
occurs for /c 11 logarthmically proportional to the inverse of the 
material loss rate, and the bandwidth in k\\ is constant. Hence, 
the radiative flux rate between the two slabs is given by: 


A ~ 47r2d2 4(Imx)^ 


(S.22) 


The asymptotic expression in Eq. (S.22) is almost identical to 
the limit in Eq. (10) in the main text, except that the flux rate 
scales logarithmically instead of linearly with |x|^/ 

Conversely, for hyperbolic metamaterials, the optimal near¬ 
field heat flux is [7] 


^(^res) 


A 
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(S.23) 


HMMs therefore do not exhibit any material enhancement; 
because the resonant modes are inside the bulk rather than at 
the surface, there is no divergence in the lossless limit. 
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III. LIMITS FOR GENERAL MEDIA 


For clarity, and with regard to practical relevance, we pre¬ 
sented in the main text only limits to heat flux between non¬ 
magnetic, isotropic bodies. Here we derive the limits for more 
general media, leading to the generalization |xp/Imx ^ 


X(Imx) X 


, as discussed in the main text. For nota- 


tional simplicity we deflne C in which case the gen¬ 


eralization is IX p / Im X 




, where 


is the 


induced matrix 2-norm [8]. This generalization applies even 
for non-reciprocal media, thanks to a generalized reciprocity 
theorem [9]. 

The Maxwell curl equations are 


V X H -j- iccD = Jg (S.24) 

— V X E + iccB = Jm (S.25) 


To simplify notation going forward, we will encapsulate elec¬ 
tric and magnetic components of flelds and currents into six- 
component vectors. We denote the flelds by t/;, the free cur¬ 
rents by cr, and the induced polarization currents by z/: 



The polarization currents within a body are related to the in¬ 
ternal flelds by the 6x6 tensor susceptibility %, 

^ (S.27) 

Given these deflnitions, the Maxwell curl equations can be 
rewritten: 


I 

1^-Vx 
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iujfiQ 


+ iujx 
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two bodies. Then the powers absorbed and extinguished by 
body 2 are given by: 


Pahs = ^ Im / v-^ 2 ^ (S.32) 

^ JVi 
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where 

^2 = -X2“' (S.34) 

Constraining Pabs < ^ext yields a limit to the absorbed 
power: 
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(5.36) 


where the second inequality follows from the deflnition of the 
induced matrix 2-norm, H-H. We can write out the squared 
magnitude of the incident held: 


IV'i 


inc,lI 


p 2 f f crt(xi)rl(x,xi)ri(x,x'i)cr(x'i) 

^ Jv^ Jv^ 


(S.37) 


The fluctuation-dissipation theorem dictates that the ensemble 
average of the current-current correlation function is 

((t(xi)(t’I'(xi)) = P [ImXi] (5(xi -Xi)0(w,Ti) (S.38) 

Inserting Eq. (S.38) into Eq. (S.37) yields the limit to the en¬ 
ergy flux into body 2 (the Planck factor separately multiplies 
the flux to give the total power): 


Following the derivation in the main text, the first step is to 
deflne a Green’s function (GF), Fi, in the presence of only 
body 1: 
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Vl JV2 


ri(xi,X2) (ImXi)r|(xi,X2) 

(S.39) 


I 

y-Vx 


Vx 

iujjj.0 




ri(x,xo) = -iujI6 (x - xo) 


(S.29) 


where it is implicit that Xi = 0 ^1 points outside of Vi. 
Then the total flelds in the presence of both bodies, excited 
by stochastic currents in body 1, satisfy the integral equation 



ri(x,xo)cr(xo) 



ri(x,Xo)X2V'(xo) 

(S.30) 


— '0inc,l 



ri(x,Xo)X2V'(xo). 


(S.31) 


Now the flelds incident from body 1 have been separated from 
the “scattered” flelds that arise only from the introduction of 
body 2, while fully accounting for interactions between the 


The integrand in Eq. (S.39) relates the flelds in V 2 , in empty 
space, from sources in Vi, within body 1. To And limits to 
this quantity, it would be useful to transpose the source and 
measurement positions in the Green’s functions. Even if body 
1 consists of a nonreciprocal material, it is possible to switch 
the source and receiver positions if the material susceptibility, 

X = I (S.40) 

VX 21 X 22 ) 

is simultaneously transformed to a complementary 
medium [9], 


Xc = 




= Sx^S 


(5.41) 

(5.42) 
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where 


5' = 



(S.43) 


and I is the 3x3 identity matrix. Defining Tc as the Green’s 
function in the presence of the complementary-medium body 
1, the modified reciprocity relation [9] dictates: 


ri(xi,X 2 ) = STq{-}C2,Xi)S (S.44) 

We can then perform a number of simplifications on the inte¬ 
grand in Eq. (S.39), including the trace operator and pulling 
the imaginary operator out front: 


IinTrri(xi,X 2 )xirl(xi,X 2 ) 

= IinTrS'r5(x2,xi)S'xiSTc(x2,xi)S' 
= IinTrS'r5(x2,Xi)xfcrc(x2,Xi)S' 

= IinTrr^(x2,Xi)xfcrc(x2,xi) 

= IinTrr^{x2,Xi)xicrc(x2,xi) 


where the first equality uses reciprocity as defined by 
Eq. (S.44), the second equality uses the definition of the 
complementary medium, Eq. (S.42), the third equality uses 
Tr SXS = Tr X, by the definition of S, and the final equality 
takes the transpose of the matrix product inside the trace. Af¬ 
ter applying these transformations, Eq. (S.39) now represents 
a new absorption problem: the absorption inside the comple¬ 
mentary version of body one due to sources in empty space 
in V 2 . This absorption problem can be bounded just as the 
previous one was, by energy conservation, such that 


(S.45) 


/ rj;.{x 2 ,xi)xic.rc(x 2 ,xi) 

Jv^ 

[ rJ(x2,Xi)ro(x2,Xi) (S.46) 

2 Jv^ 


fVi 

< 


where Fq is the free-space Green’s function and = 
—Xic' It turns out that the norm of the loss rate for the com¬ 
plementary material is equal to the norm of the loss rate of the 
original material: 


(ImCic) ^ 


-(imXic) 'II 2 

-(im [xf]“') 
- (im [Xi]”^) 


through repeated application of the facts that S ^ = S't = S' 
and that transposing a matrix does not affect its norm. Einally, 


we relate the trace of the integrand to the Erobenius norm of 
the Green’s function: 


TrTlTo = llroll^ 


(S.47) 


to ultimately yield a flux limit: 


$M<- 

TT 


(Im42) 


-1 


/ / ||ro(xi,X2)| 

' JVi JV2 

(S.48) 


that is precisely the generalization of Eq. (5) in the main text, 
for a wide class of materials. The limit could even be extended 
to inhomogeneous media, although the exact geometry would 
need to be specified to know the material loss rate everywhere. 


IV. RADIATIVE VS CONDUCTIVE HEAT-TRANSFER 
COEFFICIENTS 


We compare radiative heat transfer to conductive heat 
transfer and derive the equations used for the plots shown 
in Eig. 3(c). The total radiative heat transfer between 
two bodies is given by Eq. (1) in the main text, H = 
/ T>(co’) [©(cc, Ti) — ©(cc, T 2 )] dec. Eor a small temperature 
differential between the bodies, the conductance (heat trans¬ 
fer per unit temperature) per area A is termed the radiative 
heat transfer coefficient and is given by 

hvad = ]4 y = —ks j $(w)/(w)dw, (S.49) 

where 




hjjj 


jksT 


ksT J (^^huj/ksT _iy 


(S.50) 


When considering the limits to radiative heat transfer between 
metallic objects, one can expect that the resonances will have 
relatively small decay rates and thus that 4> will be very nar¬ 
row, and much sharper than the Boltmann-like distribution 
/(cc) in the integrand. Thus we approximate h by 


had^ ^ksfiuJo) J ^iuj)duj. 


(S.51) 


We take the metal to be a Drude metal with susceptibility 
x(cc) = + ^ 7 cc), for simplicity. Moreover, we as¬ 

sume that the absorption and emission of each body is de¬ 
scribed by a single sharp Lorentzian, with a narrow bandwidth 
(full-width at half-max) given by Acc = 7 [10, 11]. This is 
much narrower than e.g. the plane-plane and metamaterial 
structures in Eig. 3(a,b) and is in line with the resonant heat 
transfer between two spheres or between a sphere and a plate, 
depicted in Eig. 2 of the main text. The integral over 4> is then 

J <h(cc) do; = ^<h(cco) (S.52) 

and thus the radiative heat transfer coefficient is given by: 

hmd ~ (S.53) 
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The single-frequency limit to the flux per unit area is given by 
Eq. (6) in the main text, repeated here for a Drude metal: 


<h(cjo) 


< 


1 




Thus the limit to the radiative heat transfer coefficient is 

,4 


(S.54) 


ksUJo N 


(S.55) 


From a design perspective, not each of the parameters in 
Eq. (10) is a free parameter. The choice of temperature, for 
example, sets the optimal frequency (a blackbody at 300iT has 
maximum emission at 7.6/im wavelength). Similarly, the fac¬ 
tor ccp/cc is limited by the optimal aspect ratio, and the factor 
y/ccp is set by the material loss rate. Hence, it is convenient 
to rewrite Eq. (10) as 


^rad ^ 


klT 

h 


327rd^ 




(S.56) 


For a sharp, resonant spectral flux centered at cc = ccq, we can 
take 0 roughly fixed at its value at ujq, similar to the approxi¬ 
mation of / in Sec. 4. For a fiux with Lorentzian lineshape of 
bandwidth Acc, the transfer per area is 



ttAcc T>(ci;o) 


0 (cdo,T) 


(S.60) 

by Eq. (S.52). For a spectral fiux that peaks at the limit given 
in our manuscript, we have 


1 Ixl^ 

A (Imxf 


(S.61) 


where for simplicity we’ve taken Xi = X 2 = X- The Planck 
distribution factor is given by 0 = hujQl{e^ ~ 1)^ where x = 
huJolksT. For typical plasmonic systems [11] the loss rate is 
proportional to the material loss Imx/lxl^ such that we can 
approximate Acc « ujq (Imx) / |xl- Then the heat transfer 
per unit area is 


where g = j{e^ — 1)^ for x = fmo/ksT. 

The thermal conductivity of air is [12]: 

Kair = 0.0262^ (S.57) 

m • K 

Across a gap of size d, the conductive heat transfer coefficient 
is given by 


^•cond = ^ (S.58) 

d 

/irad and /icond are plotted in Fig. 3(c) in the main text for 
a variety of wavelengths and temperatures; also included are 
radiative heat transfer coefficients for plane-plane configura¬ 
tions, which fall short of the limits presented and require ex¬ 
tremely small separation distances to even reach the conduc¬ 
tive heat transfer coefficient. 


H ^ ksT Ixl" 3 a: 

A 327rc2(fcd)2 Imx^ e® — 1 

^ TT^ {ksTf Ixl^ 

4h^c^{kd)‘^ Imx — 1 


(5.62) 

(5.63) 


Two far-away black bodies exchange heat at a rate H/A = 
where ctsb is the Stefan-Boltzmann constant: 


27T^k% 

We can rewrite the near-field transfer in terms of ctsb, 


(S.64) 


f =asBT^ 


15 x^ 


ixr 


87J-3 - Ij {kdf Imx‘ 


(S.65) 


V. RESONANT HEAT TRANSFER 

One can similarly calculate the approximate frequency- 
integrated heat transfer for a narrow-band spectral fiux. The 
heat transfer is given by 

POO 

H= <I>(cd)0(cd,T) (S.59) 

Jo 

where we’ve taken one of the bodies at temperature T to be 
much hotter than the other body (such that 0i — 02 ~ 0i). 


The term in parenthesis is maximum for huj/ksT = x ^ 3.9, 
in which case the term itself is 0.28904 ... ^ 2/7, such that 
we can write 


t4 ^ 1 ^ 1 ' 


(S.66) 


which is precisely the Stefan-Boltzmann ray-optics limit, 
scaled up by the distance enhancement l/{kd)‘^ and by the 
material enhancement |xP / Ibix- 
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